Energy and multipartite entanglement in multidimensional and frustrated spin models 
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We investigate the relation between the entanglement properties of a quantum state and its energy 
for macroscopic spin models. To this aim, we develop a general method to compute energy bounds 
for states without certain forms of multipartite entanglement. Violation of these bounds implies the 
presence of these types of multipartite entanglement. As examples, we investigate the Heisenberg 
model in different dimensions, the Ising model and the XX model in the presence of a magnetic 
field. Finally, by studying the Heisenberg model on a triangular lattice, we demonstrate that our 
techniques can be applied also to frustrated systems. 
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I. INTRODUCTION 

In recent years, the investigation of entanglement in 
condensed matter systems has become one of the main 
lines of research in quantum information science Q, 0, 
3]. The increased interest in this topic is fed by several 
motivations. On the one hand, the studies helped to 
understand fundamental properties of condensed matter 
systems like quantum phase transitions. On the other 
hand, they lead to results of practical importance since 
they allowed to design new simulation techniques for the 
calculation of ground state energies of spin models 

What kinds of entanglement occur in natural situa- 
tions? This question provides another motivation for 
studying entanglement properties of condensed matter 
systems. Indeed, condensed matter systems and espe- 
cially spin models are natural candidates for our stud- 
ies where various forms of entanglement might occur, 
mainly at low temperatures. One possibility to study 
the presence of entanglement in spin systems is to re- 
late the energy or other macroscopic observables of the 
system to certain entanglement properties of the state 

In this paper we attempt to proceed in this direction by 
investigating the relation between the energy of a state 
and its multipartite entanglement properties. We will 
derive a general method for calculating energy thresh- 
olds for states without certain types of multipartite en- 
tanglement. Below these energies, and consequently be- 
low a certain temperature, the state must therefore con- 
tain multipartite entanglement. Our approach is, how- 
ever, not restricted to states in thermal equilibrium. We 
demonstrate that our method can successfully be applied 
to various models and also to frustrated systems. In this 
way, we extend the results of Ref. [10( where such en- 
ergy thresholds have been computed for two special spin 
models in one dimension. 

Our paper is divided into four sections. These are or- 
ganized as follows. In Section II we introduce the notion 



of multipartite entanglement that we use in this paper. 
That is, we explain the definition of fc-producibility. We 
also pose the problem that we want to solve. In Sec- 
tion III we present our method for computing the desired 
energy bounds. We present in detail the calculation for 
a two-dimensional Heisenberg model on a square lattice, 
the generalization to other models is then straightfor- 
ward. In Section IV we discuss three simple applica- 
tions: the Heisenberg model in various dimensions, and 
the Ising model and the XX model with a magnetic field 
in one dimension. For the Ising model, we also discuss 
the impact of phase transitions on our energy thresholds. 
In Section V we consider the Heisenberg model a two- 
dimensional triangular lattice. We show that with some 
modifications our methods can also be used to investigate 
multipartite entanglement in such a frustrated system. 



II. DEFINITIONS AND STATEMENT OF THE 
PROBLEM 



Let us first explain the notion of multipartite entan- 
glement that we use for our study. This is the so-called 
fc-producibility, introduced in Ref. |lfj |. It is defined as 
follows: For a pure state on TV qubits we ask whether 
it is possible to write 



|0l) ® |02) O •■■ ® \<j>. 



K 



(1) 



where the \<pi) are states of maximally k qubits. If this 
is the case, then only fc-qubit entanglement is necessary 
to generate \4>) and the state \ip) does not contain any 
(fc + l)-qubit entanglement. If Eq. Q holds, we call 
the state fc-producible, if not, we say that contains 
(fc + l)-partite entanglement. Examples of four- and two- 
producible states are shown in Figs. 1 and 2. 

For mixed states, we can extend this definition via con- 
sidering convex combinations. I.e., we call a mixed state 
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FIG. 1: Schematic view of a four-producible state in a spin 
model of 24 qubits. Dashed lines correspond to interactions 
between disentangled qubits and solid lines represent the in- 
teractions between qubits which are allowed to be entangled. 



g fc-producible if we can write 



(2) 



with pi > 0, pi ~ 1 and fc-producible | ipi) . If this is not 
the case, g contains (fc + l)-partite entanglement. Physi- 
cally, a fc-producible mixed state requires fc-qubit entan- 
glement and mixing for its creation only. Conversely, a 
state contains fc + 1-party entanglement if and only if the 
quantum correlations of this state cannot be explained 
by assuming fc-qubit entanglement only. 

This classification of multipartite entanglement has 
some connections to the usual notion of fc-separability, 
which is often used for small numbers of qubits. This 
issue has been discussed in Ref. ^(j- Here, we only 
want to point out that the iV-separable (fully separable) 
states are by definition the states which can be written 
as Q — Si P» Pi ® ■■■ ® 8n- These states are just the one- 
producible states. 

The notion of fc-producibility leads to a discrete classi- 
fication of multipartite states. For pure states, it is easy 
to see that the fc-producible states form a set of mea- 
sure zero in the set of 2fc-producible states and that in 
the vicinity of any fc-producible state one can find states 
with arbitrary high producibility 11]. For mixed states, 
however, this is not true anymore, and one can show (as 
for the notion of fc-separability [I^) that the set of mixed 
fc-producible states is not of measure zero in the space of 
all mixed states. 

Finally, it is worth noting that the notion of fc- 
producibility for pure states has a close relation to the 
Schmidt measure, which is an entanglement monotone 
for multi-qubit states |l3|. The Schmidt number of a 
pure multi-qubit state is defined as follows. One expands 
as the sum of tensor products of single qubit states 

IV>) = Ef=i \<Pk } ) ® \4>k ] ) ® l^ 3) ) ® ■■• For ever y quan- 
tum state we take the expansion with the minimal i?, 
which will be denoted by r. Then log 2 (r) is the Schmidt 



measure of and for N-qubit quantum states we have 
always log 2 (r 9 ) < N. 

Since the Schmidt number for two-qubit states in max- 
imally two and for three-qubit states is maximally three 
[l3|. we can conclude that for one-, two- and three- 
producible states 



log 2 ri p = 0, 



log 2 r 2p 



< 



iVlog 2 (2) N 



log 2 r 3p < ^^Uo.53iV. 



(3) 



holds. 

Now we can state the main problem we want to study 
in this paper. Let us assume that we have a macroscopic 
spin system of qubits on some lattice, which interact via 
some Hamiltonian 



H = 



E h 

<»,i> 



(4) 



which is a sum of two-qubit interactions. We always as- 
sume periodic boundary conditions. For this situation, 
we want to derive lower bounds for (H) for fc-producible 
states. That is, we want to compute a constant Ek p such 
that 



(H) > E, 



k'p 



(5) 



holds for all fc-producible states. If this bound is then vi- 
olated at low temperatures, the state under consideration 
contains (k+ l)-party entanglement. Note, however, that 
we do not restrict our attention to states in thermal equi- 
librium. Since we assume that the number of qubits N is 
large, it will be convenient to express E^ p as a rescaled 
energy per interaction bond. 

In the next Section, we will present the main idea of 
our method to compute Ek p - 



III. ESTIMATING THE ENERGY 

In this Section, we present as the main result of the 
paper a general method to estimate the energy for fc- 
producible states. While the main result is quite simple, 
its proof requires some technical effort. The logical struc- 
ture is as follows: For a pure fc-producible state some ex- 
pectation values in the Hamiltonian factorize and some 
not. We collect all the factorizing terms in the Hamil- 
tonian, and estimate them via the Cauchy-Schwarz in- 
equality. Then, we arrive at Eqs. (|13l I14fl . If we can 
perform the maximization in Eq. 113|) then the Eq. I|14|) 
delivers the desired energy bound. The point is that even 
for macroscopic fc-producible states the maximization re- 
quires only a maximization over fc-qubit states. Some- 
times, this can be done analytically, otherwise it can be 
solved numerically in an efficient manner. Finally we dis- 
cuss whether the derived bounds are sharp. 
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FIG. 2: A possible grouping for a pure two-producible state 
in twelve-qubit spin system. See text for details. 

So let us explain our method in the following example. 
We want to derive a bound for two-producible states for a 
two-dimensional Heisenberg lattice. That is, we consider 
the Hamiltonian 

Hh = ^ hij, 

<i,j> 

hij = XiXj + YiYj + ZiZj. (6) 

for a two-dimensional system with periodic boundary 
conditions. Here and in the following, Xi,Yi,Zi denote 
the Pauli matrices <7x , <?y , 0* , acting on the i-th qubit. 
We will exemplify the definitions required for our method 
using the example of a special two-producible state of 
twelve qubits given in Fig. 2. 

In order to bound the energy for two-producible states, 
it suffices to consider a generic pure two-producible state 
\ip). This comes from the fact that the mixed fc-producible 
states form a convex set in the state space and the pure 
fc-producible states are its extremal points. Thus, any 
linear function takes its maximum in a pure state as an 
extremal point. 

A fixed two-producible \ip) results in a partition of the 
whole spin system into several one- and two-qubit blocks. 
Indeed, one can identify some pairs i,j of qubits where 
the reduced state is allowed to be entangled, and some 
single qubits fc, which are not entangled with any other 
qubit. Let us denote the total number of blocks by K, 
the number of one-qubit blocks by L\ and the number 
of two-qubit blocks by L2- A possible blocking is shown 
in Fig. 2: We have seven blocks Bi,...,Bt, where B 3 
and £?5 are single-qubit blocks and the rest are two-qubit 
blocks. Thus we have K = 7, L\ = 2, and L2 = 5. It is 
important to note that we can restrict our attention to 
the case where the two-qubit blocks are between inter- 
acting qubits. This is true for the following reason: If 
a two-qubit block consists of two noninteracting qubits 
then the Hamiltonian is only sensitive to the reduced 
density matrices at each qubit, i.e., it does only take then 
local properties into account. Thus, we can replace this 
two-qubit block by two one-qubit blocks. 



In general, the mean value of the Hamiltonian con- 
sists of two-qubit expectation values of the type Wij = 
(AiAj),A = X, Y, Z. For simplicity, we denote in the fol- 
lowing (AiBj) = ciibj and (Aj) = a; for A,B = X,Y,Z. 
Note that this definition implies that in general atbj 7^ 
a, • bj and these notations have to be distinguished. Due 
to the special partition, however, some of the mean val- 
ues factorize. For instance, in the situation of Fig. 2 we 
have X2X3 = X2- X3. Now we have to find an efficient way 
for taking all these contributions into account. 

For this purpose, we define for each block Bi three sets 
of indices: are the qubits which lie inside Bi, R(i) are 
the qubits inside Bi which interact via the Hamiltonian 
with qubits outside of Bi and N{i) are the qubits outside 
Bi, which nevertheless interact with some qubit of 
For instance, in the example of Fig. 2 we have 1(1) — 
{1,2}, R(l) = {1,2} and N(l) = {3,4,5,6,9,10}. For 
our special case of a two-producible state we have = 
R(i), but in general this does not have to be the case. 

Now we can define for each block a set of expectation 
values from the Hamiltonian in the following way. We 
define 

:={W kl I k,lel{i)} (7) 

as the contributions of the Hamiltonian inside the 
block Bi. In our example, this would be 1(1) = 
\x\X2,y\y2-,z\Z2\. Note that the set T(i) is empty for 
the one-qubit blocks. In the following, we will denote 
the single elements of sets like T(%) by I(i)[j]. 

Then, we collect the "outgoing" contributions from a 
block via all the two-qubit Hamiltonians. That is, we 
define: 

K(i) := {{W kl )\ k I k G R(i),l G N(i)}. (8) 

Here, using \ k with k G R(i) we express that the two- 
qubit contributions are restricted to the operator acting 
only on the qubit belonging to R(i). For example, from 
a term of the form W23 = (X2X3) we take only the term 
(X2) = X2- Also, we take all terms with their respective 
multiplicity, i.e., if the qubit k G R(i) interacts with sev- 
eral qubits in N(i), the same term appears several times 
in TZ(i). For our example in Fig. 2 we would have: 1Z(1) = 
2/1, zi, xi, yi, z 1 ,x 1 ,y 1 ,z 1 , x 2 ,y2,z2,x 2 , 2/2, z 2 ,x 2 ,y2, 

Finally, we count for each block the contributions in 
the neighborhood via: 

:= {(W kl )\ u I k G R(i),l G N(i)}. (9) 

These are, in a certain sense, the complementary con- 
tributions to the contributions in lZ(i). We always 
write them in the same order as the contributions 
in lZ(i), i.e. the first element of TZ(i) should corre- 
spond to the first element of Af(i) in the Hamilto- 
nian, etc. In our example, we would have: Af(l) = 
{x 9l y 9 , z 9 , Xi, 2/4, Zi, x 5 , y 5 , 25, xio, 2/io, ^io, £3, 2/3, 23, x&, 
ye,ze}- The idea behind these definitions of and 
Af(i) is the following: The terms in lZ(i) and J\f(i) 
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are just the ones which factorize in the Hamiltonian. 
Thus, viewing 1Z(i) and Af(i) as vectors, the scalar prod- 
uct corresponds to the mean value of some terms in 
the Hamiltonian, TZ(i) ■ jV(i) = Efc^WM J^(i)M = 

L,k£B{i),l£B{i)(Wkl) ■ 

To estimate (H) for a given two-producible state |-0) 
we interpret TZ(i) and J\f(i) as real vectors. We then 
define 



V2 



■K(K), 



-^=-K(l)®...®-= 
' -Af(l)® ...®^=-M{K). 



V2 



V2 



(10) 



Please note that a term of the type Wm = XkXi = Xk ■ xi 
originating from an interaction between two blocks Bi 
and Bj appears twice in each of these vectors: one time 
with Xk € and xi € A/"(j) and one time with Xk € 
A/"(i) and € T^(j)- Thus, u"i and v 2 are built of the 
same terms, but in different order. This implies that 

INI = IM- 
With this definition, it follows that 



K 



i=l k 



(11) 



holds. This implies due to the Cauchy-Schwarz inequality 
that 



K 



w * EE^m-iniii^ 

8=1 fc 

- -(££-*( 



»=i /c 



(12) 



The key point is that the right hand side of this inequality 
can be estimated by minimization for each of the blocks 
Bi separately. Indeed, if we define for the block Bi 



Q :=max[^-X(z)[fc] + ^(^)[fc]) 2 ], (13) 

k 



2 ^ 

k 



where 



is a quantum state on the block Bi, we have 
(H)>-J2d- (14) 



The estimation of the Cj does now only depend on the 
fact whether the block Bi is a one- or a two-qubit block 
and not on the relations between these blocks. For the 
Heisenberg interaction, we have 



Ci = max [2(xl + yl + zf)] 
m 



(15) 



for a one-qubit block Bi on the qubit k and 
I 

+ \{xl+yl + zl + xi+yf + zf) 



max [ - x k xi - y k yi - z k zi + 
\i>) 



13 



(16) 



for a two-qubit block on the qubits k and /. This bound 
can be obtained from the representation of = 
Sfc i=i x y z ^ki&k ® &i ,1Q|. It is also a special case of 
a general bound presented as Lemma 1 in the Appendix. 

With these bounds, we immediately get for our exam- 
ple in Fig. 2 the bound (H H ) > -(5-13/3+2-2) = -77/3. 
For the general case of N qubits, we get 



,„ v r „ r 13 T , UN 
(H H )>- max 2-Li + — -L 2 > — 

Li+2L 2 =AT L 3 6 



-2.16iV. 
(17) 

Since a two-dimensional lattice of N qubits has 2N 
bonds, the energy per bond for two-producible states is 
bounded from below by 



(Hh) 
2N 



> E 2p — — - 



13 
12' 



(18) 



Two questions arise at this point. First, we have to 
ask whether this bound is useful, in the sense that it is 
violated at low temperatures. This is the case since for 
the ground state the energy per bond is Eq = —1.338 
Thus, in a considerable temperature regime, the 
thermal states cannot be two-producible. 

Second, the question arises whether the derived bound 
is sharp. This question deserves some discussion. The 
idea to show sharpness of an obtained bound is the fol- 
lowing: Let us assume we have found states \<fii) for 
which the maxima Ci in Eq. I|13l) are obtained. Then 
we have to build out of these states \4>i) the total state 
\4>) such that \ip) saturates the bound in Eq. (jnj. To 
do so, we have to assure, that for the state \ip) the 
Cauchy-Schwarz inequality in Eq. (|12|) was sharp, i.e. 
Vi-V2 — — ||^i ||||^2 1|- This can be done in two steps: First 
we guarantee that v\ -v% = \\v\ \\ \\v2W ■ Then, by applying 
some unitary transformations on the \<f>i) we make sure 
that vi ■ v 2 = -||ui||||# 2 ||- 

Let us show how this works in our example of two- 
producible states 0|- Let \<j>)k,l be the state satu- 
rating Eq. i|13fl on the qubits k,l. Let us enumerate 
the qubits as in Fig. 2. and consider the total state 
\i>) = |0)i,2 ® W)a,4 <8> |0)5,6 ® |0')7,8 <» |0)9,io-- Here, 
|0') = 5(|0)) is the state |</>) where the qubits are 
swapped. This construction implies that the reduced 
states of on the qubits 1,4,5,8... are identical, as 
well as the reduced states on the qubits 2, 3, 6, 7... Since 
the corresponding reduced states are identical, the fac- 
torizing terms between two qubits (say, 2 and 3) are just 
squares of some expectation values, hence V\ and v 2 are 
parallel and v\ ■ v 2 = \\vi\\ 

To perform the second step, note that the state |0)fc,z 
on the qubits fc, I gives rise to some sign distribution of 
the expectation values Xk,yk,Zk and xi,yi,zi. Then we 
define \<j>") as follows. We first swap, i.e. |0")m,n = 
|0')m,n — <S(|0)m,n)> then, by local unitary transforma- 
tions, we flip the signs of x m , z m and x n , z n on the qubits 
m, n. Finally, we transpose the density matrix of the 
state, which flips also the signs of y m and y n . Thus we 
have finally cik = —a n and ai — —a rn for a = x,y, z and 
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Note that \<f>") still saturates Eq. i|16|) . since the ex- 
pectation values a m a n are not affected. Then, defining 
\ip) = \<f>)i,2 <& |0")3,4 ® |0)5,6 ® ■•■ we arrive at a state 
for which v\-V2 — — ||wi||||w2||. Thus, this state saturates 
Eq. JH|). 

In general, however, the bounds derived by the method 
above are not sharp. Especially, when we consider frus- 
trated lattices, the bounds are not sharp, and more so- 
phisticated estimates are required. We will discuss one 
example of a frustrated lattice later in detail. Also, if N 
is not a multiple of fc, the bound for fc-producibility may 
not be sharp. This is, however, not a major problem. 
The energy difference between this case and the nearest 
N which is multiple of k is bounded by a constant. Since 
we are interested in the thermodynamic limit N — ► 00, 
and the energy difference per bond decreases as 1 /N, we 
can neglect this case. 

What is required to derive similar bounds as Eq. i|18|) 
for other spin systems and higher degrees of multipar- 
tite entanglement? The main ingredient are bounds as 
in Eqs. 1151 IT^f) . These bounds depend on the Hamil- 
tonian and on the underlying lattice. For many Hamil- 
tonians and two-qubit blocks, these bounds can straight- 
forwardly be computed analytically. But even if this is 
not possible, one can simply compute them by numerical 
minimization over a small number of qubits, if desired, 
this minimizations can be performed with assurance of 
global optimality [l6[ . 

Finally, the reader should note the difference between 
the estimation method presented in this Section and the 
one used in Ref. 0. The method in Ref. does not 
separate between factorizing and non-factorizing terms 
in the Hamiltonian, instead, it estimates the complete 
Hamiltonian via the Cauchy-Schwarz inequality. As a 
consequence, it requires more effort and is restricted to 
one-dimensional systems. 



IV. THREE SIMPLE APPLICATIONS 

In this Section, we will apply the presented method to 
various examples of spin systems. We will first compute 
energy bounds for fc-producibility of spin systems with 
an anti-ferromagnetic Heisenbcrg interaction in various 
dimensions. Then we will consider the Ising model and 
the XX model in an external magnetic field. 



Then, the energy bounds per bond for one- 
and four-producible states are given by 



two-, three- 



V 1D 



1: 



^3p 



-1.505; 



77.113 

&2p 
&4p 



3 

-1.616. 



(20) 



The ground state energy per bond is known to be Eq = 
-(4 In 2 - 1) « -1.773 17], thus all the bounds above 
are violated by the ground state. 

(b) For the two-dimensional square lattice, the respective 
energies per bond are given by 



El? = -1; 
E™ = -1.108; 



T? 2D — 

&2p - - 

W 2D _ 



13 
12' 
-1.168. 



(21) 



Here, the energy per bond in the ground state is given 
by E = -1.338 0. 

(c) For the three-dimensional lattice we have Eq — 
— 1.194 [3 and the thresholds for multipartite entangle- 
ment read 



E 



3D 



- ^2p - ~ — ' 



£§f = -1.044; E\» 



31 
30' 
-1.067. 



(22) 



All the bounds given in this theorem are sharp. 

Proof. The proof of this theorem works just as de- 
scribed in the previous section. For the one-producible 
(i.e. the fully separable) states, the bounds have al- 
ready been shown before 0, IE • The bounds for two- 
producible states have been obtained analytically (see 
Lemma 1 in the Appendix); here, the bound for the one- 
dimensional chain was already derived in Ref. jlCj . The 
bounds for the three- and four-qubit case have been ob- 
tained numerically ^8|- Note that for four-producible 
states and D > 2 several possibilities of four-qubit blocks 
have to be taken into account. The sharpness of the 
bounds follows also as discussed in the previous section. 
□ 



B. The Ising model in a transverse magnetic field 

As a second example, let us study the one-dimensional 
Ising-model in a transverse magnetic field. That is, we 
consider the Hamiltonian 



A. The Heisenberg model 

Let us first consider the anti-ferromagnetic Heisenberg 
interaction. That is, we consider the Hamiltonian 



H 



H 



E 

<ij> 



•V , A , + V, V, + /., /., 



(19) 



on a D-dimensional lattice. For this model, we can state: 
Theorem 1. (a) Let us consider an infinite one- 
dimensional spin system with the Heisenberg interaction. 



H 



(23) 



<h0> 



The estimation of the energy for fc-producible states can 
be performed as in the previous section. Only the inter- 
action terms with the magnetic field have to be added in 
the definition of X(i). For instance, for two-qubit blocks, 
we have to compute 



C, 



1 



max 



x k xi - B{z k + zi) + - (x k + 



(24) 



0.5 1 1.5 2 2.5 3 3.5 4 
B 

FIG. 3: Entanglement gap E g (k, B) for two- three-, four- and 
five-partite entanglement for the Ising model in a transverse 
magnetic field. See text for details. 

This and similar maximizations can easily be performed 
numerically. The resulting bounds are always sharp. 
Note that the bound for one-producible states has al- 
ready been derived employing a different method in 
Refs. H0. 

The Ising model is analytically solvable and the ther- 
modynamic properties of the thermal states are known 
[l9|. To investigate the multipartite entanglement prop- 
erties, we first compute the energy thresholds Ek p for k- 
producible states. We then compare these energies with 
the ground state energy by calculating the entanglement 
gap 

E g (k,B)=E kp (B)-E (B), (25) 

that is the difference between the ground state energy 
and the minimal energy for fc-producible states Q ■ Note 
that the energy minimum for separable states for a quan- 
tum Hamiltonian equals the energy minimum of the cor- 
responding classical spin chain jfj. Thus E g (l,B) is the 
energy difference between the classical and the quantum 
Ising spin chains. The results are shown in Fig. 3. 

To discuss these results, let us consider Fig.[3]and look 
at the curve corresponding to E g (l,B), For a magnetic 
field slightly larger than B c = 1 the entanglement gap 
(and thus the entanglement in the thermal state) is larger 
than further from this point. Note that at B c — 1 the 
ground state of the Ising model undergoes a quantum 
phase transition. E g (n, B) for n > 1 also takes its maxi- 
mum around B c . Fig.[3]shows that the field corresponding 
to this maximum is decreasing with increasing n. 

Now let us study F(k,B) = dE g (k,B)/dB as the 
derivative of E g (k,B) with respect to B. These curves 
are shown in Fig. 21 On can see directly from this fig- 
ure that the slope of F(1,B) has an abrupt change at 
B c i = 2. Further analysis shows, that F(1,B) is also 
non-analytical at B c = 1. The non- analytical point B c 
corresponds to a quantum phase transition of the quan- 
tum spin chain, while the change in the slope at B c \ = 2 



-0.1 
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B 

FIG. 4: Derivative of the entanglement gap F(k, B) = 
dE g (k, B)/dB for k = 1, 2 and 4. 



corresponds to the critical point of the classical spin sys- 
tem |20j. 

Similar results can be obtained for E g (k, B) for k > 1. 
The energy minimum for fc-producible states equals also 
the energy minimum of a spin model in which blocks of 
k quantum spins interact classically, i.e., in a mean- field 
fashion p]| . Again there are non-analytical points for 
F(k,B) at B c = 1 and at B c ^ n > 1. Fig. 01 shows the 
curves corresponding to F(1,B), F(2,B) and F(4,B). 
It is clearly visible how B c k approaches B c = 1 with 
increasing n. It can also be seen that the maximum of 
F(k, B) also approaches B c as k increases. A detailed 
study of the thermodynamics arising from these models 
intermediate between classical and quantum spin chains 
will be reported elsewhere. 



C. The XX model in a magnetic field 

As a third example, we study the one-dimensional XX 
model in a magnetic field. The Hamiltonian of this model 
is 

Hxx = X * X i +Y 1 Y. j +bY j Z 1 . (26) 

<i,j> i 

The estimation of the energy for fc-producible states can 
be performed similarly as for the Isi ng m odel. Since the 
XX-model can be solved analytically |23| , it is now inter- 
esting to investigate the regions in the T-S-plane where 
multipartite entanglement must be present. This has 
been done in Fig. 5. Similar to the Ising model after the 
quantum phase transition at B = 2 the thermal states 
show different forms of multipartite entanglement, even 
at relatively high temperatures. 
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not 4-producible 




not 2-producible 




not 3-producible 




not 1 -producible 





FIG. 5: Entanglement in thermal states of the XX-model 
in a magnetic field. The regions in the T-B-plane are shown 
where the different types of multipartite entanglement can be 
detected with our method. 



V. THE HEISENBERG MODEL ON A 
TRIANGULAR LATTICE 



FIG. 6: A two-producible state on a triangular lattice. Solid 
lines represent possible entanglement between the qubits. The 
triangle (2, 5, 3) is type A, i.e., it does not have entanglement 
between its qubits. The triangle (1,2,3) is of type B since 
qubits 1 and 2 may be entangled. When estimating d for 
the block of the qubits 1 and 2, the two triangles (1, 2, 3) and 
(1,2,4) are estimated via Eq. (13 XP . since there are definitely 
of type B. See text for further details. 



Let us finally demonstrate with an example that our 
method with some modifications also allows the compu- 
tation of energy thresholds for frustrated lattices |23| . 
Generally, all lattices can be divided into two classes: 
Bipartite lattices are lattices, where the lattice points 
can be divided into two sublattices, such that each point 
in each sublattice interacts only with points which be- 
long to the other sublattice. An example is the two- 
dimensional square lattice, for which these two lattices 
form a chessboard-like configuration. A lattice is called 
frustrated if it is not bipartite. This terminology refers to 
the fact that for such lattices the ground state energy per 
bond is usually larger than that for two qubits interacting 
alone. 

Entanglement properties of frustrated systems have 
also been investigated 0, 13 • Concerning our approach, 
the fact that the ground state energy is large makes it 
difficult to derive energy bounds for fc-producible states 
which are violated by the ground state. 

As such an example of a frustrated quantum system we 
study now the Heisenberg model on a two-dimensional 
triangular lattice. That is, we consider the Hamiltonian 
of Eq. JBJ on the lattice of Fig. 6. Let us shortly note 
some properties of this system. The ground state energy 
per bond in known to be Eq = —0.726 [L4J. From a com- 
parison with a classical spin configuration, it was shown 
in Ref. that the minimal energy per bond for fully sep- 
arable (i.e. one-producible) states is E\ v = —0.5. Here, 
we want to derive a bound for two-producible states. 

If we apply directly the method of the previous section, 
the resulting bound is not violated by the ground state. 
The reason is the following: In the derivation, we used in 
Eq. I|12fl the bound v\ ■ V2 > — \\vi\\ \\v2\\- This bound is 



not sharp for frustrated lattices. Thus, we have to make 
a more sophisticated estimate. 

First, note that the scalar product v\ ■ V2 represents 
all factorizing terms in the Hamiltonian. These terms 
can be grouped into the contributions corresponding to 
different triangles T{. So we can write 



2V X -V2= 

triangles T; 



(27) 



where the triangle contributions rj(Xi) can be of two 
types, depending on the triangle: For a triangle on the 
qubits j,k,l with no entanglement between the qubits 
j, k, I (type A triangle, see Fig. 6) we have 



t)(Ti) = ^ (aj ■ a-k + a k ■ ai + ai ■ aj). 



(28) 



For triangles where two of the three qubits (say, k and /) 
may be entangled (type B) we have 



h(Tj) = ^ (a, • a k + dj ■ ai). 



(29) 



The prefactor of two in Eq. l|27Jl stems from the fact that 
every bond contributes to two triangles. 
Now we need the facts that 

/ a,j ctj a k ak ai ai \ 

(30) 

which holds for all real numbers aj , ak, fflj, an d we need 
the estimate 



'akai 



1 a,- 



a k a t 



1 a 



y/2 V2 y/2 



< ak ■ aj +ai-aj. (31) 
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This estimate holds since and a k ai are expectation 
values of (tensor products) of Pauli matrices. Namely, 
from the positivity of the density matrix the relation 
— (1 + afc<2;) < afc + ai < (1 + aicO-i) follows, which re- 
sults in Eq. (3TJ|. 

The bounds in Eq. Ij30l I31|l can be interpreted in the 
following way: for each of the f)(Ti) there are two vectors 



TO and u>„- 2 ^ such that 



-(1) -(2) 
-W\ ■ W\ . 

-(1) 



(32) 



If Tj is of the type A then w\ has nine entries and if Ti 
is of the type B it has six entries. Note that the definition 
of vfp etc. does not directly correspond to the 1Z and 
Af in the definition of Vi/ 2 . 
Thus, if we define 



tC 1 ) 



(33) 



it follows that 2v x ■ v 2 > ~W {1) ■ W (2) . Using the Cauchy- 
Schwarz inequality yields 

vi-v 2 >-^(\\W^\\ 2 + \\W^\\ 2 ). (34) 

The right hand side of Eq. (|34|l is a sum of many 
squares of expectation values of one- or two qubit observ- 
ables. Each of these expectation values originates from 
a certain block. By counting carefully the contributions 
of each block, we can now estimate them separately for 
each block. 

A one-qubit block on the qubit fc, contributes to the es- 
timates of six triangles. These may be triangles of type A 
or B. An estimate of one triangle of the type A results in 
two times a contribution (a fe /V2) 2 in ||JU (1 )|| 2 + ||JU (2) || 2 , 
one in ||1UW|| 2 and one in ||IU( 2 )|| 2 (see Eq. (JSOJ). A 
triangle of the type B gives also two times this contribu- 
tion, but now either two times in ||W'^ 1 ' > || 2 or two times 
in ||VK( 2 )|| 2 (see Eq. JSTJ). Thus, in analogy to Eqs. (EH 
I15L we have to estimate 



(35) 



A two-qubit block on the qubits k and I contributes 
to ten triangles. Two of them contain both the qubits k 
and I and are thus of the type B (see Fig. 6). Each of 
them contribute two times (a k ai + l/-\/2) 2 to HJU^H 2 + 
||IU (2) || 2 , either both in \\W^\\ 2 or both in ||IU^|| 2 (sec 
Eq. (|31[n . For the other eight triangles, it does not mat- 
ter as for the one-qubit blocks whether they are of the 
type A or B. Thus, we have to estimate 



Ci = max — XkXi 

2-2Al + x k x l f 
4 ' 2 



VkVi ~ z k zi + 



{1 + VkVif , {\ + z k zif 



4 1 2 



vl 



2 



This follows similarly as before, see also Lemma 2(b) in 
the Appendix. 

With these bounds we have for the two-producible 
state (H) > — . Ci , which finally results in a minimal 
energy per bond of 



E- 



2p 



(37) 



(36) 



which is clearly violated by the ground state. Thus, at 
a considerable temperature, in a frustrated triangular 
Heisenberg lattice the thermal state is not two-producible 
and the effects of multipartite entanglement cannot be 
neglected. 

Finally, it should be noted that it is not clear whether 
the bound in Eq. (|3Tjl is sharp. Two facts suggest that 
this is not the case. First, the use of the Cauchy-Schwarz 
inequality in Eq. (|34|l is probably not sharp, since all 
the vectors and are usually not parallel. Sec- 

ond, the maximum in Eq. (|36|l is obtained for a separable 
state, no entanglement is required to saturate this bound. 



VI. CONCLUSION 

In conclusion, we have developed a method to inves- 
tigate the presence of multipartite entanglement in spin 
models. Our method relies on energy bounds for cer- 
tain classes of multipartite entangled states and can be 
applied to any state, not only to thermal states. We 
discussed different examples and showed that our ideas 
can in principle also be applied to the investigation of 
frustrated systems. 

The results of this paper may be applied in several di- 
rections. On the one hand, based on our energy thresh- 
olds one may derive temperature bounds below which 
entanglement must be present. Then, one may try to de- 
sign methods to extract this entanglement and make it 
useful for some tasks. 

On the other hand, our results can also be used to gain 
theoretical insight concerning the validity of ground state 
approximations. For example, one might be interested in 
the ground state energy Eq. To approximate this, one 
may consider some mean-field like approximation, where 
the trial wave function is a product state with respect 
to all qubits [2!|. This corresponds to a one-producible 
state. Concerning the justification of this approximation, 
we can say with our method the following: If E\ p « E 2p 
the approximation might be justified, while if E\ v E 2p 
this is clearly not the case. 
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VIII. APPENDIX 



yed j g bright to diagonal form via a singular value de- 
composition. This decomposition can be performed by 
local unitary operations. 

Thus it suffices to consider \if>) = a\00) + /3|11) with 
a 2 + [3 2 — 1, since for that state we have 



M = 



1 a 2 -(3 2 

2a(3 

-2a(3 
a 2 - (3 2 1 



(40) 



Here, we prove two useful bounds for our proofs. 
Lemma 1. For two-qubit state on the qubits k,l 
and 7 > 1/2 we have the following sharp bound: 

-x k XL - y k yi - z k zi + j(x 2 k + y\ + z\ + xf + yf + zf) 

<l + 2 7 +^. (38) 
27 

Proof. We have to maximize the left hand side of the 
inequality over all states. First note that a generic quan- 
tum state can be written as 



\ijO~i 



(39) 



with o~i — 1. Here, Xij — {ip\<Ji ® o~j\ip) holds, thus we 
can directly maximize over all allowed Xij. Let us first 
consider local unitary transformations. These transfor- 
mations act on as (Ay) — ► (1 ffi Oi)(Ajj)(l © Or) 
where Ol and Or are orthogonal 3x3 matrices. Here, 
1 © Or denotes a 4 x 4 matrix with a block structure, 
i.e., with "1" in the left upper corner and Or denotes the 
3x3 block in the right bottom corner. These transfor- 
mations do not change the purities of the reduced states, 
thus they do not change x\ + yf + zf and xf + yf + zf. 
Furthermore, IccfcX;! + \y k yi\ + \z k zi\ is the sum of the 
absolute values of the diagonal elements of the 3x3 sub- 



matrix AJ 



Xi j—x 



This sum is maximized when 



The final maximization over all a can then be directly 
performed. The left hand side of Eq. (|38l) is maximized 



for 4a 2 = 2 - 
Eq. 



I/7 which proves the bound in 

□ 



Lemma 2. (a) For two-qubit state \ip) on the 
qubits k, I and 7 > 1/2 we have the following sharp 
bound: 



-x k xi-y k yi+j(x] 
(b) Similarly, we have 



+ 7{x k +y k + xf + yf) < l + 2 7 + 



1 

87' 



(41) 



\ + \i x kxf + y k yf + z k zf) + 

+ {x 2 +y 2 +z 2 +x 2 +y 2 +z 2 )<A. 



(42) 



Proof, (a) For the minimization it does not matter 
whether we minimize over the observables X, Y [as in 
Eq. I)4ip] or X, Z. Then, the bound can be derived as in 
the proof of Lemma 1. One arrives again at Eq. 14011 . 
now only one of the terms 2a/3 has to be omitted in the 
final maximization, (b) This can also be derived as in 
Lemma 1. When applying local unitary transformations, 
((x k xi) 2 + (y k yi) 2 + (z k zi) 2 ) is maximal, when the 3x3 
matrix is diagonal. Finally, one has only to maximize 
over a again. □ 
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